Topological field theory of dynamical systems 
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Here, it is shown that the path-integral representation of any stochastic or deterministic 
continuous-time dynamical model is a cohomological or Witten-type topological field theory, i.e., 
a model with global topological supersymmetry (Q-symmetry). As many other supersymmetries, 
Q-symmetry must be perturbatively stable due to what is generically known as non-renormalization 
theorems. As a result, all (equilibrium) dynamical models are divided into three major categories: 
Markovian models with unbroken Q-symmetry, chaotic models with Q-symmetry spontaneously bro- 
ken on the mean-field level by, e.g., fractal invariant sets (e.g., strange attractors), and intermittent 
or self-organized critical (SOC) models with Q-symmetry dynamically broken by the condensation 
of instanton-antiinstanton configurations (earthquakes, avalanches etc.) SOC is a full-dimensional 
phase separating chaos and Markovian dynamics. In the deterministic limit, however, antiinstan- 
tons disappear and SOC collapses into the "edge of chaos". Goldstone theorem stands behind 
spatio-temporal self-similarity of Q-broken phases known under such names as algebraic statistics 
of avalanches, 1/f noise, sensitivity to initial conditions etc. Other fundamental differences of Q- 
broken phases is that they can be effectively viewed as quantum dynamics and that they must 
also have time-reversal symmetry spontaneously broken. Q-symmetry breaking in non-equilibrium 
situations (quenches, Barkhausen effect, etc) is also briefly discussed. 

PACS numbers: 05.45.-a, 05.65.+b, 02.50.Fz 



Many fundamental aspects of the theory of non- 
linear dynamical systems have already been es- 
tablished. Nevertheless, our understanding of 
nonlinear dynamics may benefit from turning to 
alternative approaches. Here, it is shown that 
one of such alternatives is Witten-type or co- 
homological topological field theories (W-TFTs). 
fll-llOl] It turns out that any dynamical model in 
its path-integral representation is a W-TFT. This 
approach allows to clarify the physical essence of 
Intermittency also known as self-organized criti- 
cality (SOC, for a review see, e.g., Refs.[llJ) and 
establish its connections with the other two fun- 
damental concepts of Chaos and Markovianity. 



I. INTRODUCTION 

For the last two decades, mathematical physicists and 
mathematicians have formulated and have been devel- 
oping mathematical constructions known as topological 
field theories (TFTs, for a review see, e.g., Ref.flf). TFTs 
come in two types: Schwartz-type or quantum TFTs and 
Witten-type or cohomological TFTs (W-TFTs). Quan- 
tum TFTs already found their applications - they are 
believed to stand behind exotic low-temperature phases 
of condensed matter systems such as fractional quan tum 
Hall effects and superconductors (see, e.g., Ref.pjJ and 
Refs. therein). As to the W-TFTs, to the best of our 
knowledge they are still used only for purely mathemat- 
ical purposes. From the discussion in this paper it will 
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follow that W-TFTs is actually the path-integral version 
of the dynamical systems theory. The later, in turn, has 
many applications in modern science. 

It is well known that the most general Parisi-Sourlas- 
Wu (PSW) stochastic quantization procedure [13|, [l4| 
applied to Langevin equations leads to N=2 (quasi- 
)Hermitian supersymmetric models (Witten models, see, 
e.g., Refs. [l|, [l5|). It is also known that determin- 
istic conservative (classical, Hamilton) dynamical sys- 
tems viewed as path-integrals also possess supersymme- 
try. fig, [l7[ The above supersymmetries are of topo- 
logical origin. [J] In this paper, it is demonstrated that 
topological supersymmetry is pertinent to any determin- 
istic or stochastic continuous-time dynamical model in its 
path-integral representation. Based on this and on the 
possibility of the spontaneous breakdown of the topo- 
logical supersymmetry by two different mechanisms, a 
generic phase diagram for dynamical models is proposed. 
The phase diagram is given in FigHjD. 

As quantum field theories, W-TFTs that show up on 
PSW quantization are non-Hermitian. Therefore, our 
proposition relies strongly on recent developments in the 
theory of non-Hermitian quantum dynamics. [18l420l | 

The paper is organized as follows. In Sec|Hl it is 
demonstrated that any stochastic or deterministic con- 
tinuous dynamical model is a W-TFT. In Sec lIIIi we fo- 
cus on models with white noise. Non-Hermitianity of the 
models (Sec lIIID|) . their spectrum and physical states 
(Sec lIIIE|) . etc., are addressed. In Sec lIVl the meaning 
of the ground states in deterministic limit is analyzed. It 
is shown that in cases of unbroken topological symmetry 
the collection of (bra's) ket's of the perturbative ground 
states is the representation of the (anti-)instantonic CW- 
complex of the phase space. In turn, the global ground 
gma^tatfls represent stable and unstable manifolds that inter- 



sect on invariant manifolds. Sec. [V]is devoted to models 
with unbroken topological symmetry, and, in particular, 
to Langevin models. It is also discussed how W-TFT's 
with no physical quantum excitations represent stochasti- 
cally fluctuating dynamical models (Sec|VB]). In Sec lVIl 
models with spontaneously broken topological symmetry 
are addressed. Qualitative difference of dynamics with 
spontaneously broken topological symmetry is discussed 
in general terms in Sec lVIAl Particularly, it is proposed 
that such concepts of dynamical systems theory and com- 
plexity theory as self-similarity, sensitivity to initial con- 
ditions, non-Markovianity etc. are actually mere con- 
sequences of the Goldstone theorem. Two mechanisms 
of spontaneous topological symmetry breaking are iden- 
tified: on the mean- field level, e.g., in the determinis- 
tic limit, and due to the dynamical condensation of in- 
stantons and antiinstantons. It is argued that these two 
mechanisms corresponds respectively to chaos fSec lVIBJ) 
and intermittency/SOC fScc lyTC]) . In Sec. IVTDl non- 
equilibrium dynamics such as quenches or Barkhausen 
effect is briefly addressed in the context of W-TFTs. Sec. 
IVIII concludes the paper. 



II. STOCHASTIC QUANTIZATION 

Dynamical systems are those that are defined by spec- 
ifying their equations of motion. This is the most general 
class of models as compared to, say, the equations of mo- 
tions that follow from least action principle. 

The equations of motion can be either step-like equa- 
tions or time-continuous equations. This paper deals 
with time-continuous evolutions only. Continuous dy- 
namics can be referred to as physical dynamics since time 
is always continuous in physical systems. 

The equations of motion can be either stochastic or 
deterministic (partial) differential equations (SDE or 
DDE). DDE's can be looked upon as SDE's with zero- 
variance noises. Therefore, it suffices to construct a path- 
integral formulation for an SDE, and the theory of a 
corresponding deterministic model will follow by sending 
the parameters of the noise to their deterministic limit. 
Thus, SDEs are of primary interest. 

The procedure of building a partition function out of 
an SDE is called stochastic quantization. There are two 
maj or stochastic quantization procedures: PSW method 
PS EH and Martin-Siggia-Rose (MSR) method. [H| As 
is discussed below, MSR is an approximation to the most 
general PSW method that we use in this paper. 

Consider an SDE in its most general form: 



F i (R,<p)=?(R,<p). 



(1) 



Here R s B with B being some Riemannian base space 
having the meaning of spacetime and ip S M(E, B), with 
M being an infinite-dimensional space of all maps from B 
to a topological target manifold, S. In Eq. ([1]), F's have 
the meaning of a DDE and £'s is the stochastic noise. F's 



are some functionals of </?'s and may have explicit depen- 
dence on the base. Both F's and £'s belong to the tangent 
space of M, i.e., F l (R,tp)6/5tp l (R),C(R,tp)6/5tp l (R) G 
T^itmM. F l, s and £'s can also be thought of as sections 
of the tangent bundle TM . 

The noise must obey an important physical condition. 
It must experience no feedback from the system. If there 
is such a feedback, the noise should be viewed as a part of 
the system rather than as an external source of stochas- 
ticity. This condition, however, does not necessarily sug- 
gest that the stochastic correlators of the noise have no 
functional dependence on ip. For example, below we con- 
centrate on models where the noise is decoupled from the 
system but its correlators are functions of (p. This pos- 
sible dependence of the correlators of the noise on ip's is 
emphasized by the explicit dependence of the noise on tp 
in the rhs of Eq.(JTJ. 

The above condition only suggests that there exists an 
invertible transformation of TM — > TM: 



(<P,Z) -> (<P,£), 



(2) 



that decouples the stochastic variables from (p's. An ex- 
ample of such a transformation is given in Eq. (|16[) . The 
partition function of new stochastic variables is: 



((1))- // 1-e 



-s„ 



=(«) 



(3) 



where S no i ae is independent of ip's. 

In Eq. ((31) and in the following the double integral sign 
denotes functional integration. ((...)) is a common nota- 
tion denoting a vacuum expectation value in the theory 
of quantum non-Hermitian models, which the models we 
consider here will turn out to be. We accept this notation 
from the very beginning of our discussion. 

Now, the original SDE ([1]) can be rewritten as: 



F a (R,&=£ a (R), 



(4) 



where F a (R, (p) = £, a (<p, F(R, tp)) is the functional trans- 
formation from Eq.@ applied to F's. Eq. (fTT)) below is 
an example of the transformed SDE (0| . 

We assume that the noise is "decent" and for any, 
X a {R), there exist a unique and well-defined functional 

S- \ sc {X) = \o g {(e^^ x ^)) 

oo „ n 

= £/ C^ ) - a "(R 1 ...R n )l[X ai (R i )/n\, (5) 

where C's are the irreducible stochastic correlators for 
£'s. Here and in the following the factors representing 
the metric on B are assumed whenever appropriate. 

The inverse functional transformation in Eq. ([5]) is 
well defined, e.g., when the noise has only one "vacuum" , 
i.e., only one solution for 5S no i se /5t; a (R) = 0. This vac- 
uum must be stable and non-degenerate so that all the 
eigenvalues of S 2 S no i SC /6^ a (R)5^ b (R') are positive. Such 



" instanton-free" noises can be called Markovian in the 
context of this paper. 

The hrst step toward a W-TFT, which to the best of 
our knowledge belongs to Parisi and Sourlas, [13j is very 
natural. It is the realization of the fact that the only 
partition function a stochastic system may in principle 
have is that of its noise, Eq.©. Indeed, the partition 
function is the summation over all the realizations of a 
stochastic process, which is actually £. Furthermore, if 
we want that the " number" of the noise's degrees of free- 
dom be the same as the number of the system's degrees 
of freedom we must consider only temporally periodic 
boundary conditions. This is nothing else but the equi- 
librium situation that we will consider mostly. 

The second step is to rewrite the partition function in 
terms of (p's rather than £'s. Let us now recall that in 
the theory of supersymmetric models there is a concept of 
Nicolai maps. [l|, [22J It says that for theories with global 
supersymmetry one can come up with such a transfor- 
mation of bosonic variables, that the partition function 
will fold into that of a noise because the Jacobian of the 
variable's transformation will cancel the fermionic deter- 
minant. What we need to do now is essentially the oppo- 
site. We have to unfold the partition function of the noise 
into some other model and rightfully anticipate that the 
model will possess a global supersymmetry - topological 
supersymmetry (Q-symmetry). 

The partition function of the noise is rewritten as 



provides the desired Jacobian 



«1» = 



m] 



5{F a -C)Je 



-S„ oJ se(f) 



(6) 



where J = det(SF a (R, tp)/S<p b (R')) is the Jacobian of the 
variable transformation. The (5-functional limits the in- 
tegration over M only to the solutions of the SDE, that 
serves as an "inverse" Nicolai map. In order to bring 
Eq.© back to Eq.©, one must integrate out £'s and 
notice that the outcome is actually Eq.© with £'s sub- 
stituted by <p's according to the variable transformation 
specified by Eq. © . 

The Jacobian of the variable transformation is impor- 
tant. Without it, Eq.© is not the same as Eq.© and 
thus does not represent the stochastic process under con- 
sideration. Neglecting the Jacobian is the mathematical 
essence of the MSR stochastic quantization procedure. 
Therefore, MSR picture can only serve as an approxima- 
tion, while the PSW stochastic quantization procedure, 
the one we use here, is always correct. 

In the standard fashion, one introduces the Lagrange 
multiplier, B a , to incorporate the (5-functional: 



S(F a -£ a ) 



„i! R B a {F«-t) 



(7a) 



Here and in the next formula, ~ denotes equality up to 
an unimportant constant. One also introduces the set of 
the Fadeev-Popov ghosts, x a Xb, integration over which 



J 



,-if Ra >X*(SF a /6 l p b ) x b 



[\x] 



(7b) 



On inserting Eqs.© into Eq.© one obtains 



«1))=// e^/^C^-nl-^tt). (8) 

J ■/[£*] 

Here and in the following $ stands for the collection 
of all the fields CpBxX- These fields constitute what 
could be called the supersymmetric extension of TM. 
In Eq. © , we introduced the operator of the gauge-fixing 
Becchi-Rouet-Stora-Tyutin supersymmetry, which is the 
Q-symmetry: 



{Q,x} 



X a (R)S^ {R) +B a (R)S Xa{R) \X. (9) 



Here and in the following 8 (or d) denotes functional 
(or partial) differentiation over its subscript. Operator 
Q is a bi-graded differentiation, {Q,XY} = {Q,X}Y + 
(— l) m X{Q, Y}, where m is the ghost degree of X. Q is 
nilpotent, {Q, {Q,X}} = 0, for any X. 

Out integration of the noise in Eq. © leads to 



«1» = 



J{QJ r x*f"}+s- o \ so {-i{Q,x*}) 



(10) 



where the last term in the action is defined in Eq.©. 
This term is a functional of only Q-exact pieces (of the 
form {Q,X}) and is Q-exact itself because of the nilpo- 
tencyofQ: {Q, X}{Q,Y}... = {Q,X{Q, Y}...}. Thus, 



((1)) = 



*{S,e} 



(11) 



with the so-called gauge fermion given by 

e = / xa(R)F a (R) + J2(-*T +1 / XaARi) 

JR n=2 JRx...R n 

n 

xC^ ) ~ a "(Ri...R n )l[B ai (R i )/nL (12) 



The last term in this expression comes from the noise 
and for deterministic models it vanishes. 

A Q-exact action is the unique feature of W-TFTs that 
all look like a gauge fixing of "nothing". That the PSW 
method leads to a W-TFT comes with no surprise. The 
point is that we are computing the partition function of 
the " ignorant" external noise, which does not really care 
about the lhs of Eq.©. Therefore, the deformation of 
the SDE must not result in any changes in the partition 
function - a feature pertinent to W-TFTs. 

Topological nature of models with Q-exact actions and 
with non-trivial topological content can be in particu- 
lar revealed by the provided possibility to calculate on 
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FIG. 1: (a) Graphical explanation of how negative probabil- 
ities show up in stochastic quantization. Original noise has a 
positive delta-function-like probability distribution, (b) Cor- 
responding to a highly nonlinear many-to-one variable trans- 
formation, F(if) = £ (dashed curve), that has the mean- 
ing of an SDE, the probability distribution, P(tp) (solid 
think curve), is negative in those regions where the Jacobian 
dF/difi < 0. 



instantons certain topological invariants as expectation 
values of collections of certain Q-closed operators also 
known as BPS (Bogomornyi-Prasad-Sommerfield) oper- 
ators. 10J Nevertheless, even W-TFTs with topologically 
trivial content are of topological origin - in Ref. [4| , any 
model with a Q-exact action has been given an interpre- 
tation of a generalized Morse theory on M. Some new 
mathematical insights on W-TFTs have been recently 
developed in Ref.j 



A. Connection to Negative Probabilities 

Stochastic quantization seems to have a close connec- 
tion with the concept of negative probabilities (see, e.g., 
Ref.[23| and Refs. therein). This concept shows up nat- 
urally on the PSW quantization. 

Consider, e.g., a nearly-deterministic Gaussian noise, 
^gK 1 , ina (0+0) theory. The partition function is: 



((1)) 



P(0 = / P(v) = * 



= *{Q,e} 



(13) 



Here, P(£) ~ e - ^ / 2e with e — > 0, the second integral 
is over a new variable that is related to £ through some 
many-to-one function, £ = F((p), P((p) = P(F((p))J(<p) 
is the new probability density with J(ip) — dF(<p)/dip 
being the Jacobian of the variable transformation, and 
the last equality sign establishes the connection to the 
discussion in the previous subsection with © — x(F(ip) + 
ieB/2). 

New probability distribution is given in Fig.([T|). In 
those regions where the Jacobian is negative, P(<p) is 
also negative. These negative probabilities are necessary 
for the partition function to be equal to the partition 
function of the original stochastic process. They always 
appear paired with new positive probabilities and thus 
are clearly of topological origin - the simplest realization 
of the Poicare-Hoft theorem. 

A mathematician would probably say that P(ip) is a 
pullback of P(£) by the irreversible many-to-one map, 



F(ip) : ip — ► £. The original P(£) is non-negative and 
thus can be interpreted as a volume form for £. From 
this point of view, the stochastic quantization is the pro- 
cedure of "borrowing" the volume from the noise, so that 
the total volume in the new variables (the total proba- 
bility to exist) remains the same as in Eq. (|13p . 

This simple example suggests the following stochas- 
tic interpretation of negative probabilities. If a stochas- 
tic variable has negative probabilities, it means that its 
partition function actually represents/calculates a parti- 
tion function of yet another "physical" stochastic vari- 
able (noise) with ordinary positive probability. The two 
stochastic variables must be related through some irre- 
versible many-to-one map, e.g., through a highly nonlin- 
ear SDE. 

On the side of stochastic dynamics, our temporary un- 
derstanding of the negative probabilities is this. For a 
specific configuration of the noise, a nonlinear SDE may 
have a multitude of solutions that in the high energy 
Physics terms are called Gribov copies. The system has 
a freedom to chose which one of the Gribov copies is go- 
ing to be realized. We can not know which one of the 
Gribov copies the stochastic system is going to chose. 
Only this chosen copy must contribute to the partition 
function. At the same time, we have to integrate over 
all the possible configurations of the systems variables 
(over the entire M) and all the Gribov copies are going 
to contribute. This would necessary lead to overcount- 
ing unless some of the Gribov copies are contributing -1 
instead of 1 due to that the probability density is nega- 
tive in the corresponding regions of M. Thus, negative 
probabilities represent the freedom of a stochastic system 
to chose among various competing solutions of its SDE. 
From the forthcoming discussion, it follows that this free- 
dom gets physically realized only in Q-broken phases. 



III. WHITE NOISE CASE 

For the following discussion we do not need the gener- 
ality of the previous section. We assume that spacetime, 
R = (t, x) 6 B = T x S , where T is time and the space, 
S, is flat (e.g., a torus S = T d , i.e., periodic boundary 
conditions). In equilibrium situations T — S , while in 
non-equilibrium cases T = R . We discuss mostly the 
equilibrium situations. 

SDE has only the simplest temporal non-locality: 



w+a^c 



(14) 



where the flow, A 1 , is some functional of tp at this specific 
moment of time. A can be though of as a vector field over 
the (infinite-dimensional) phase space, H, of all maps 
from the space to the target: H = {(p % (x) : S — » £}. 
The noise is assumed Gaussian and white: 

((e(R)?(R'))) = 5{R-R')g l \ V {R)), (15) 

where g lJ has the meaning of the metric on E, which is 
a function of only tp 1 's at this specific R. 



Eq. (fT4|) and/or Eq.flll]) with Eq. (|19b[) below represent 
a very wide class of stochastic dynamical models. In 
particular, this class includes systems defined by SDE's 
that are not first-order in time-derivative. One can bring 
such higher-order SDE's to the form of Eq. (fT4)) by the 
introduction of new fields. In case of Kramer' s equation 
this works like this: dfip — A = £ is rewritten as dtp~A = 
£, dt^p—p = C , where £' is a zero- variance Gaussian noise. 

The transformation that decouples the noise from </?'s 
has the following form 



r(i?) = e n<p(R))c(R), 



(16) 



where e's have the meaning of vielbeins: e°5 a f,e^ = cjij 
g %0 . The transformed SDE from Eq. d3| 



and e l J ah e\ 

takes the following form: 



e<}{d t ^ + A l )=t- 



(17) 



Straightforward application of the PSW stochastic quan- 
tization procedure from the previous section to the par- 
tition function of the white noise: 



«1» 



,-I R it?/2 



[?] 



(18) 



leads to a W-TFT for the collection of fields, $ = 
(ip\x\B a ,Xa), with a Q-exact action, S = {Q, 6}, de- 
fined by the gauge fermion, = J R Xa(ef(d t ip' 1 + A 1 ) + 
iS ab Bi,/2), and the supersymmetry operator {Q,X} — 
J R (X 1 (R)5^ {R) +B ? (R)5^ {R) )X. 

It is more convenient, however, to introduce a new set 
of fields, $ = (ip l ,x l ,Bi,x.i), where ip l and x % are the 
same, while Xi = X a ef and B t = {Q,Xi} = B a ef - 
Xk^uX 1 with rf ; = e„(e")J being the Christoffel symbol. 
There are two things to note in this redefinition of the 
fields. First, {Q,B t } = {Q, {Q,Xi\} = automatically 
because of the nilpotency of Q. Second, the measure in 
the pathintegral in unchanged. 

In terms of the new fields, the model is a W-TFT 



<0» = // e l{S ' e} , 



(19a) 



with 



6 = / Xl {prf + A 1 + igVtBj + x*lV)/2) , (19b) 
Jr 



{S,X}= / (x l (R)6^ iR) +B l (R)S MR) )X. (19c) 
Jr 

We purposely used the transformed SDE method to 
quantize the model in order to make direct connection 
with the discussion in the previous section. In fact, model 
(fT9| can be obtained with less effort (with no redefinition 
of the fields) - out integration of the noise in the model 
constructed from the original SDE (fT4|) : 



«1» 



e HQ,f R X,(3t</-K4 s -eir)}-(f) 2 /2 



[«*] 



(20) 



with Q from Eq.(|19c|). leads directly to Eqs.([I9|). 

For d = 0, model (|T9| is topological quantum mechan- 
ics, [4| in which case H = S. In many cases, below we 
consider this situation. Higher dimensional theories with 
d > can be viewed as infinite-dimensional generaliza- 
tions of topological quantum mechanics. 



A. Schrodinger picture 

In passing from the path-integral representation of the 
theory to the operator algebra representation, the parti- 
tion function is rewritten as: 

{{i))=[[ e f t (L(*B z d tV *-ixid tX i-Hm)^ ( 21 ) 

From here, the composition laws for the operators follow 
in the standard manner 

[^(x),4-(:r/)]- = -\x i (x)M*')]+=ti i A*-*'n 22 ) 

The subscripts denote commutation for bosonic fields 
and anticommutation for the ghosts. Other (anti- 
commutators are zero. 

We chose to work in the representation where ip's and 
x's are diagonal. There are two reasons for this choice. 
First, these fields are superparthners. Second, in this 
basis the Noether charge associated with the topological 
symmetry is the exterior derivative that has no explicit 
dependence on the metric, which emphasizes the topo- 
logical nature of the model. 

In this basis, Eq. (|22"T) suggests 



Bi(x) = -iS v i( x ),x.i(x) = ~i8 x i( x ). 



(23) 



Wave functions, |\1/), are functionals of ip's and x's only. 
We can formally Taylor expand in x's: 

„=ol*> (n) ' 
l*) (?l) = ^)...( Xnkn )*(x k Hxi)...X K (xn)).(2A) 

Here, the integration over x's and summation over fc's 

(n) 

is assumed, and w^ , -. , , ■> are antisymmetric in the 
pairs of subscripts due to the anticommutation composi- 
tion law for the ghosts. For this reason, |\I/)( n )'s can be 
interpreted Q as forms from the exterior algebra of H 
which thus is the Hilbert space, H, of the model. 

In Eq. ([24|) . * denotes the Hedge star, which is intro- 
duced for convenience. Within this definition, the prob- 
ability density, which is a form of maximal degree, cor- 
responds to |*)(°) as in Eq. (j2?|) below. 

The dynamical equation governing the time-evolution 
follows directly from the path-integral formulation in 
Eq. ([2~T|) and is the generalized Fokker-Plank equation: 



d t m 



-H\V). 



(25a) 



Its time reversed version for bra's, (\1/| = (*|^ r ))*, that 
are also forms on H, is 



d t (V\ = {^lijt. 



(25b) 



The Hamitlonian in Eqs. (j25)) is the generalized Fokker- 
Plank Hamiltonian. It's explicit form can be derived 
straightforwardly by the bi-graded symmetrization of its 
path- integral expression and with the use of Eq. (l2"5)l . The 
result is well-documented in the Literature [10( | and is 
know to have a form of a N=2 (pseudo-)supersymmetric 
(pseudo-)Hermitian model: [19| 



H =[Q,Q}+/2 = -\7 2 /2-C A . 



Here 



Q = X l {x)8 vHx)l 



(26a) 



(26b) 



is the conserved Noether charge associated with Q- 
symmetry. Operator Q is nothing else but the exterior 
derivative on W. Besides Q, the Hamiltonian also con- 
serves the number of ghosts given by the operator: 



In Eq.flHaD, 



F= / x'O^x'C*)- 



t -2 A\x)5 xi{x) , 



(26c) 



(26d) 



is what could be called the (pseudo-)conjugate super- 
charge with 



& 



^■WS^^W-xr^^,)), (26e) 



being the adjoint of the exterior derivative. 

In Eq. (|26ap , — V 2 = [Q, Q^]+ is the Laplace-Beltrami 
operator, explicit form of which is provided by the 
Weitzenbock formula, and Ca = [Q, f A*(x)5 x ir x \]+ is 
the Lie derivative along A. 

We would also like to note, that formally the model 
is not a N=2 pseudo-supersymmetric in general for that 
reason that Q is not nilpotent, Q 2 ^ 0, and it is not 
commutative with the Hamiltonian, [H, Q]_ =/= 0. 



B. Connection to conventional Fokker-Plank 
equation 

Because number of ghosts is the integral of motion, the 
time evolution does not mix wavcfunctions of different 
ghost number. This, in particular, suggests that the time 
evolution of the wavefunction of the trivial (or rather 
maximal) ghost content: 



depends only on itself: 

dtP = -H cnv P, (28) 

where (g l ' 2, s come from the Hedge star in Eq. (j2"7| '): 



rrenv 



~ 1/2 W9 1/2 (V J 'W)/ 2 - Ai ( x )) >( 29 ) 



is the conventional Fokker-Plank operator acting on the 
probability density, P. 

If from a physical point of view, the probability density 
is a function(al), P, from the topological point of view, 
it is a form of the maximal degree as it is in Eg. (1271) . 
Notably, in stochastic quantization the probability den- 
sity is not a " square" of a wave- function (of trivial ghost 
content) but rather the wave-function itself. The cor- 
responding bra is actually the cycle over the entire H. 
This unconventional meaning of wave- functions is the di- 
rect consequence of the non-Hermitianity of the model 
that we discuss below 



C. Meaning of wavefunctions and Q-symmetry 

In relation to the physical meaning of P, it is worth 
to briefly discuss the physical meaning of wavefunctions 
with non-trivial ghost content. In literature, |\t)W is 
interpreted as current. 24 [4J| It represents the flow 
of probabilities. This interpretation can be supported 
within the TFT formalism in the following manner. Con- 
sider topological quantum mechanics with d-dimensional 
phase space. Consider also Stock's theorem: 



Jen Jn 



1^(^,1) =^P(^,t)), 



(27) 



where f2 is a (fixed, time-independent) d-dimcnsional 
part of the phase space, <9f2 is its (d — l)-dimensional 
boundary, and l^}' 1 ' is any wavefunction with d — 1 
ghosts. We know that (Q\^)^) has the meaning of prob- 
ability density so that \^)^ must have the meaning of 
the probability flow though d£l. Inductively, wavefunc- 
tions of even more non-trivial ghost content must have 
the meaning of currents of currents of probability etc. 
Within this interpretation of wavefunctions, the model 
complies with the continuous higher-dimensional version 
of Kirchhoff's law. [25| Then, operator from Eq.(26d) 
must be identified as conductance (or rather as conduc- 
tance over capacitance). 

The previous Stock's equality must hold at any mo- 
ment of time and thus it shows that I^W and Q]^)^ 
evolve in time equivalently. Time evolution for both is 
given by Fokker-Plank equation so that the Fokker-Plank 
Hamiltonian must be commutative with Q, which is the 
Nother charge of Q-symmetry. In other words, if a the- 
ory describes a dynamical model in terms of probability, 
currents of probabilities etc. it must possess topologi- 
cal supersymmetry in order to be consistent with Stock's 
theorem. 



This meaning of wavefunctions also suggests that sta- 
tistical description (in terms of probability density) of 
stochastic models is not applicable when Q-symmetry 
is spontaneously broken. Indeed, Q\ ground) ^ auto- 
matically suggests that the ground state is not a (pure) 
probability density since <3|*) (0) =0 for any |<J} (0) . 

Wavefunctions with non-trivial ghost content may also 
be given the alternative interpretation of conditional 
probability densities. Even though there are a few phys- 
ical reasonings that support this interpretation, we will 
not pursue this line of thinking in this paper. 



D. Non-Hermitianity 

There is a natural composition for bra's and ket's: 



<*l|*2> 



(**l)* A * 2 



(30) 



H 



In Hermitian/unitary models, Eg. (1301) is the metric in 
the Hilbert space, %. In models under consideration, 
however, Hamiltonians are not Hermitian and the metric 
is not the one provided by Eq. (I31H) . To see how non-trivial 
metric shows up in non-Hcrmitian models, we address 
here the spectrum and eigenstates of the model. 

It is convenient to think of H as of an infinite dimen- 
sional matrix with real entries. Its spectrum consists 
of either real eigenvalues or pairs of complex conjugate 
eigenvalues. If complex conjugate pairs exist, the model 
is said to be pseudo-Hermitian, while if the spectrum is 
purely real the model is said to be quasi-Hermitian. [l9( 

For pseudo-Hermitian Hamitlonians, there exist a bi- 
orthogonal basis in the Hilbert space: 



H\n}} = £ n \n)),((n\H = ((n\£ n . 



(31) 



Here, \n)) = \n) are merely the eigenstates of H, while 
((n\ = (m\fj mn are related to (n|'s through some non- 
trivial metric of the Hilbert space such that, 

fjHfj- 1 = #+. (32) 

The bi-orthogonal basis is complete: 

((n\m)) = 6 nm , i n = Y] \n))((n\, (33) 

where 1% is the unity operator on H. The metric mixes 
the "original" eigenstates with complex conjugate eigen- 
values because: 



(n\fjH = (n\fjS* 



(34) 



The pairs of states with complex conjugate eigenvalues 
are time-reversal companions. This means that time- 
reversal operation acts non-trivially on such states. 
The dynamical equation for a generic ((^>\ = (&\fj is 



instead of Eq.(|25bp. 

In the literature on W-TFTs, pathintegrals as the one 
in Eq. (|2T|) . are often called Euclidian. This identification 
has a hidden danger in the context of dynamical systems. 
It may sound like the " actual" (or physical) temporal evo- 
lution is governed not by Eq. ([25l) but by the Schrodinger 
equation related to Eq.(|25|) through Wick rotation of 
time. The point here is that the time in Eq. (f2"Tj) is the 
actual time of the SDE and the fundamental dynami- 
cal equation is the Fokker-Plank equation as it should. 
Therefore, apart from the subtle question of the validity 
of Wick rotation in pseudo-Hermitian models, this line 
of thinking may lead to an accidental "switching" of the 
physical meanings of the real and imaginary parts of the 
eigenvalues of the Hamiltonian. 

To separate the real and imaginary parts of the eigen- 
values we introduce: 



&n — *■ r, 



•iE„ 



(36) 



dt((*\ = {{*\H, 



(35) 



In stochastic quantization, T's are the attenuation rates 
(inverse lifetimes of eigestates), while E's are quantum 
mechanical energies (inverse periods of closed classical 
trajectories associated with the eigenstates). This in par- 
ticular suggests that stochastic quantization has a subtle 
difference with quantum dynamics. The role of the con- 
ventional kinetic term in quantum systems (~ —A/2) is 
to set the lower limit on the energies thus ensuring the ex- 
istence of a ground state. In stochastic quantization, the 
kinetic term sets the lower limit on the attenuation rates 
thus caring more about the "stability" of the dynamics. 

From physical reasonings it follows that the lower limit 
on T's must be zero. Indeed, the states with T < 
would grow infinitely as time flows, thus indicating some 
sort of instability in the pathintegral representation of 
an SDE. On the other hand, the pathintegral represents 
actually the stochastic noise, which does not have any 
instabilities. Thus, states with T n < must not exist. 
On the other hand, the partition function (which is again 
that of the noise) can not vanish in the T — > oo limit 
and thus must possess at least one "physical" state (see 
below) with T being exactly zero. 

The above physical reasoning about the possible 
generic form of the spectrum has its counterpart in the 
dynamical systems theory (see, e.g., Sec. Ill of Ref . [27J ) . 
There, the trace of the evolution operator plays the simi- 
lar role as the partition function in the W-TFT approach. 
Perron- Frobenius operator is the analogue of the Fokker- 
Plank operator of the W-TFT. The zero-eigenvalue state, 
which must always exists, is called "ergodic" zero. All the 
other eigenstates must obey T > 0. Complex eigenval- 
ues also come in conjugate pairs and are called Ruelle- 
Pollicott resonances. 



E. Physical and ground states 

Let us now turn back to the path-integral represen- 
tation of the theory. Consider the expectation value of 



some observable, O: 



«0» 



0e i{C,0}_ 



[<I>] 



(37) 



We are considering the equilibrium case, so that all the 
fields have periodic boundary conditions on a large tem- 
poral circle t £ [T/2, —T/2]. By integrating out the fields 
Si's and x's one obtains Eq. ([37|) in the operator algebra 
representation: 

((0))=Tr n (-lfT[O s U], 

where Os is the observable in Schrodinger picture, U — 



-J t i 



-T/2 



H(t) 



is the time evolution operator, and T de- 



notes chronological ordering. In the previous formula, 
(— 1) F shows up because of the periodic boundary con- 
ditions for ghosts (see, e.g., Ref. J3l|). In "normal" case 
of anti-periodic boundary conditions for fermionic fields, 
this factor does not exist. 

One can now plug the unity from Eq. (1331) at temporal 
infinity, ±T/2, and rewrite Eq. (|3"T|) as: 

{{O)) =J2(-l) Fn ((n\T{d S U}\n)), 

n 

The operator of temporal evolution acts on bra's and 



ket's in a simple manner: U\n)) = e 



-£ n T 



n)),((n\U 



((n\e £nT . Using this, we further rewrite Eq. (|3"7l) as: 



((o)) = E(- 1 ) 



F„ -£„T 



((n\T[6 H }\n)}. 



(38) 



Here, T[Ojj] = U T\&sU\ is the chronologically or- 
dered O in the Heisenberg representation. Diagonal ma- 
trix elements of T[Oh] do not depend on T. For ex- 
ample, if the operator is a two-time correlator, O = 
X(t)Y{t'), t > t', its diagonal matrix element is 



((n\T[6 H }\n)) 



-(£ m -£ n ){t-t') 



(39) 



Now it follows that in the T->0 limit, only states with 
T rl = contribute to Eq. (|37l) . Only these states can 
be considered physical ones as only they survive the in- 
finitely long Fokker-Plank evolution and thus can appear 
as the out-states in the scattering matrix. 

In fact, r„ = should be looked upon only as a neces- 
sary condition for a state to be a physical. Other condi- 
tions may apply. One of the examples of such conditions 
exists in quantum field theories where states that have 
negative norm are viewed as non-physical. This condi- 
tion does not apply, however, to W-TFTs (see, e.g., Sec. 
3.6.3. of Ref.[l|). Our temporal understanding of the 
physical meaning of the negative-norm states (those with 
odd number of ghosts) is through the concept of negative 
probabilities discussed in Sec. Ill Al 

We can adopt yet another condition for physicality of 
a state by requiring that the norm of a state is non-zero. 
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FIG. 2: (a) Graphical representation of the manipulations 
with the path-integral in Sec. (|IHE|) . Curly arrows represent 
the direction of time. The infinitely far future and past are 
integrated out around some finite time domain t > t' , of 
some observable denoted as a collection of dots. This evolves 
the eigenstates according to ((n\(t) — ((n\e~ n ^ T ' 2 ~ t ' and 



!»»(*') 



-£„(t'-(-T/2)) 



n)). In the physical limit, T — > oo, 



only the physical states with zero real parts of their eigenval- 
ues, En,, survive, (b) Qualitative comparison of the general- 
ized Fokker-Plank Hamiltonian spectra in models with unbro- 
ken (top) and spontaneously broken (bottom) Q-symmetry. 
E n — lm£ n and F n = Re£ n are the imaginary and the real 
parts of the eigenvalues of the Fokker-Plank Hamiltonian, £ n . 
Unphysical states with F n > 0, physical states with T„ — 0, 
and physical ground states with minimal E n are represented 
respectively as hollow circles, circles with grey filling, and 
black circles. Q-symmetry is broken on states with £ n ^ 
and time-reversal symmetry must be broken on states with 
E n 7^ (Ruelle-Pollicott resonances). In case of unbroken 
Q-symmetry, the only physical states are the ground states. 



We will encounter such possibility in the next subsection. 
At this moment, we are not aware of any other conditions 
for states to be considered physical. For this reason, allow 
us to call all states with T„ = physical ones. The 
Hilbert space H is split now as: 



U = W®U U , 



(40) 



where W is spanned by all the physical states, and "H u 
is spanned by the non-physical states with r„ > 0. 
In the T->oo limit, Eq. (|38|) takes the form 



«°»= E (- 1 )' 



iE n T 



((n\T[6 H }\n)). (41) 



|n»SW 



One notes that there is an oscillating factor, e~ lEnT , that 
suppresses the summation (or rather the integration in 
the higher dimensional theories) over the physical states. 
In the spirit of quantum field theories, one can formally 
Wick rotate "a little" the parameter T -> T(l - *0 + ). 
Again, in the T — > oo limit this will eliminate all the 
physical states from Eq. (|4T|) except those with the lowest 
Eg = min„ e -HP E n . These physical states could be called 
ground states. In result we get: 



«0»~ E (-^((n\T[d H }\n)). 



(42) 



\n))eH<> 



where T-L 9 is the part of the Hilbert space spanned by the 
ground states. 

The meaning of the similarity sign in the last expres- 
sion is twofold. First, we dropped the factor e~ %E a J ', 
The second meaning is more fundamental. The point is 
that for Q-symmetry broken cases, Eq. (|42|) may not be 
a good approximation for some observables. An example 
of such an observable is the unity operator. For O = 1-n, 
the average is the Witten index 



«i}} = Ec- 1 : 



F ne -T£ n 



(43) 



which is a well known fact that the partition function of 
a W-TFT is the Witten index. The only difference in the 
stochastic quantization case is that £ n 's can be complex. 

As it will be discussed in the next subsection, all the 
eigenstates with non-zero eigenvalues break Q-symmetry 
and thus must come in bosonic-fermionic pairs, who's 
contribution must cancel each other out of the partition 
function. Therefore, only the zero-eigenvalue states must 
contribute. In result, in topological quantum mechanics 
Witten index equals Euler characteristic of the target 
manifold, which in general is not zero. This must be true 
even for cases with broken Q-symmetry - after all, the 
partition function is that of the noise and thus it must 
not depend on the flow. On the other hand, if we used 
Eq. (j42|) for the calculation of the partition function in 
this cases, we would get zero. Consequently, it is Eq. (|4"T) 
that must be used for a general observable. 

Turning back to the path-integral representation 
brings Eq. ([4Tj) to a "more" field-theoretic form: 



«°» = E (-!) J 

|n»ew 



{{n\Oe 



{S,e} 



n}}, (44) 



that will be helpful in the next subsection. 

The purpose of the textbook-level exercise in this sub- 
section was twofold. First, it showed that the non-trivial 
metric of the Hilbert space of a pseudo-Hermitian model 
is automatically incorporated into the pathintegral rep- 
resentation. [26| Second, it revealed the essence of the 
unphysical, physical, and ground states of the model. 



F. Q- and time reversal symmetries of eigenstates 

Q-symmetry is pertinent to path-integral representa- 
tions of all dynamical systems. A generic eigenstate, 
however, does not possess the symmetry of the Hamilto- 
nian. If this is true for a ground state, it is said that the 
symmetry is broken spontaneously. Therefore, in order 
to judge if Q-symmetry is broken or not, it is important 
to understand which eigenstates break it. This is the 
primary goal of this subsection. 

The most general criterion of that an eigenstate does 
not break some continuous symmetry is that the expec- 
tation value of any observable, which is a result of acting 



by an infinitesimal operator of this symmetry on some- 
thing, is zero. For Q-symmetry, this criterion takes the 
following form (c.f., Eq.(j44j): 



(45) 



((n\{Q,X}e i ^ & y\n))=0, 



for any X. Using integration by parts with respect to Q 
and {Q, e^ s ' e >} = 0, Eq.® reduces to: 



{Q, |n))} = Q\n)) = 0, 
±{Q,({n\} = ((n\Q = 0. 



(46a) 
(46b) 



Here we used the fact that bra's and ket's depend only on 
ip, x so that Q operator from Eq.Q becomes Q operator 
from Eq. (|26bj) . The choice of sign in the second line 
depends on the ghost degree of ((n\. 

In combination with Eqs. (f3"Tj) and (126a[) . Eqs. (l46l) also 
suggests that: 



QQ\n)) = 2£ n \n)), 
{{n\QQ = 2£ n ((n\ 



(47a) 
(47b) 



with Q from Eq. (|26dp . Accordingly, bra's and 
ket's of eigenstates that do not break Q-symmetry 
and at the same time have £ n ^ are Q-exact: 
\n)) = Q\x n )),((n\ = ((x n \Q, where \x„)) = 
Q\n))/{2£ n ),{{x n \ = {(n\Q/{2£ n ). This, in turn, sug- 
gests that {{n\n)) = ((x n \Q 2 \x n )) = 0. In other words, 
if Q-symmetry is not broken by an eigenstate with non- 
zero eigenvalue, then it has zero norm. This situation is 
exotic. If such situation occurs, some measures must be 
undertaken to get rid of this state. 

The temporary conclusion is that all the eigenstates 
with non-zero eigenvalue break Q-symmetry. From this it 
follows that the Q-symmetry is unbroken spontaneously 
if the only physical states (those with r„ = 0) are those 
with E n — 0. In other words, the only physical states are 
the ground states that satisfy (n 6 H 9 = H p ) 



Q\n))=Q,((n\Q = 0, 
QQ\n)) = 0,((n\QQ = Q. 



(48a) 
(48b) 



These conditions are very similar to those of the unbro- 
ken N=2 pseudo-supersymmetry [19J with the only differ- 
ence that the second line is weaker than its N=2 pseudo- 
supersymmetric version: Q\n)) = 0, ((n\Q = 0. 

The comparison of spectra of models with broken and 
unbroken Q-symmetry is given in FigOa. In essence, the 
spontaneous Q-symmetry breaking can be interpreted as 
the condensation of (a branch of) Ruelle-Pollicot reso- 
nances (i.e., eigenstates with imaginary eigenvalues) into 
the physical part of the Hilbert space. 

Speaking of time-reversal symmetry, as it follows from 
Eq. (|34|) . each Ruelle-Pollicott resonance, either con- 
densed or not, is not the time-reversal companion of it- 
self. This suggests that time-reversal symmetry must be 
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broken on these states. It then follows (see Fig|2}j) that 
the spontaneous breakdown of Q-symmetry must always 
be accompanied by the spontaneous breakdown of time- 
reversal symmetry. This line of thinking is very similar to 
that about the spontaneous breakdown of N=2 pseudo- 
supersymmetry in Ref.|20], where it was shown that it 
must always come together with the spontaneous break- 
ing of time-reversal symmetry. 



IV. DETERMINISTIC LIMIT: GROUND 
STATES 

In this section, we would like to discuss the form 
of the ground states of topological quantum mechan- 
ics ((0+1) theory) in the deterministic limit for cases 
of unbroken Q-symmetry. The discussion will provide 
a route to understanding the origin of the spontaneous 
Q-symmetry breaking on the mean-field level, which hap- 
pens in chaotic models (see Sec lVIBp . 

In the deterministic limit, when the Fokker-Plank 
Hamiltonian reduces to the Lie derivative, the ground 
states, |0)), are closed states (Q|0)) = 0) that satisfy 



C A \0)) = 0. 



(49) 



Lie derivative, however, is not a bounded operator. 
Therefore, a more appropriate route is to perform a one- 
loop analysis of a model with the metric: 



<r 



in 



(50) 



and then send the 



ministic limit, e — > 0. The exact form of Qn must not 



'intensity" of the noise to its deter- 

matter as long as it is decent, e.g., positive definite. 

Let us first consider an isolated critical point, 
(pO,A % (ipo) = 0. The one-loop partition function of the 
tower of the one- loop eigenstates around ipo is: 



«|l|»L l0 ° P 



A&Mvo)} 






F ne -T£ n 



(-l) A -o. (51) 



Here the gauge fermion generating the corresponding 

gaussian action is 0(<£>o) = Jt=_T/2 A'i^t'V + AjSipi + 

ieg^Bj), Sip 1 — tp l — ip l , n enumerates the one- loop eigen- 
states, A) = d V jA l (ip ) so that A l (tp) w A) 5<p>. 

The last equality sign in Eq. ([5~Tj) follows from the fa- 
mous cancelation of bosonic and fermionic determinants 
in W-TFT's (see e.g., Ref.Q), with A Vo being the num- 
ber of negative real eigenvalues of A 1 , . 

Eq. (|5"Tj) says that there exist (at least one) one-loop 
state with zero eigenvalue of the Hamiltonian. The zero- 
energy state is either bosonic or fermionic depending on 
A^. All the states with non-zero eigenvalues come in 
bosonic-fermionic pairs so that their contribution cancels 
out from the partition function. 



In this subsection, we call zero-energy state a ground 
state. As long as we are talking about the unbroken Q- 
symmetry case, this is reasonable. 

We consider what could be called the minimal "non- 
potential" generalization of a potential flow near tpQ. 
We assume that there exist local orthogonal coordinates, 
ip l = O'jSipi , in which A 1 , is block-diagonal with diago- 
nal elements being either 1 x 1 or 2 x 2 matrices. One- 
dimensional blocks corresponds to real eigenvalues of A 1 , 
and to the directions, in which the flow is locally poten- 
tial. The two-dimensional blocks correspond to pairs of 
complex conjugate eigenvalues that represent, e.g., sinks. 
All eigenvalues have non-zero real part. Otherwise, the 
critical point is not isolated and belongs to a higher- 
dimensional invariant manifold. 

In these coordinates, the directions corresponding to 
different eigenvalues of A^ become independent (we chose 

g l Q that does not mix coordinates of different diagonal 
blocks of Aj ) . The ground state will be the wedge prod- 
uct of the ground states in each of the coordinates or the 
pairs of coordinates: 



-i-r(A) -i-r(A) 



L /9 



(52) 



where a and j3 numerate respectively the coordinates of 
the single real eigenvalues of A 1 , and the pairs of coor- 
dinates of the complex eigenvalues. Similar factorization 
holds for bra of the ground state. 

We consider first a coordinate, Cp a = ip, corresponding 
to a real eigenvalue of A 1 -, A. We chose metric g^ * = 
1. It is straightforward to verify that the ground state 
satisfying Eqs. (|48l) with: 



Q = xd^,Q = d x 



-ed^ - 2\<p), 



is 



l"» 



-\\\v 2 /t 



xAH 



(53) 



(54) 



for a stable coordinate (A > 0) and \a}} O ((a| for an 
unstable coordinate (A < 0). 

In the Hermitian basis of a Langevin model (see 
Eqs. (|63b[) and (|63cp below) and with the corresponding 
superpotential, W = \<p 2 /2e, the wave-function localizes 
to the critical point 



\Ot L 



*(Hl)* 



\\\0 2 /(^) 



X, A>0, 
1, A<0, 



(55) 



This is the well-known one-loop ground state in Witten 
models (see, e.g., Chap. 10.3.3 of Ref.[31(). The only 
difference here is that we are using the representation 
in which x's are diagonal instead of x's, so that x o 1 
brings Eq. (|55p to its more conventional form. 

Let us turn to two coordinates, (p° = {(p l ,(p 2 ), with 
a complex conjugate pair of eigenvalues, A = A' ± A". 
First, we note that the two coordinates contribute the 
factor +1 into Eq. flBTI) . This suggests that the ground 
state is bosonic, i.e., it has either two or zero ghosts. As 
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in case of a single eigenvalue of A, the first case corre- 
sponds to two stable coordinates (A' > 0) and the second 
case corresponds to two unstable coordinates (A' < 0). 
Let us consider the first case. 
Without loss of generality: 



A\ 



A' 



A"e" 
A' 



(56) 



with k > 0. We can now use the freedom in Eq. (|5T)|) of 
choosing a metric and set 



,9o 



i+/A,A=(A-A+) 1 / 2 . 



(57) 



Here A+ = (i + A T )/2 and A ± = A' ± A" sinh k are the 
eigenvalues of A + . At this, the condition shows up that 
A + is positive definite. We assume it is satisfied. 
It can be straightforwardly verified that: 

l/3)}~e- A ^/Vx 2 ,«/?l~l, (58) 

with ip 2 = (if 1 ) 2 + {(p 2 ) 2 , satisfies Eqs.gS]) with: 



Q = x l £>^,Q = d x ,( 



-£9od^ 



2A*0>) 



(59) 



where i, j — 1,2, and go and A from Eqs. (l57|) and (|56|) . 
In a similar way, for a pair of unstable coordinates one 
similarly gets the ground state with |/3)) -f-> ((f3\. 

Eqs. (fMf and (j5"5)) show that in the deterministic limit, 
e — >• 0, and under certain conditions, the one-loop bra 
and ket of the ground state of an isolated critical point 
localize and/or represent respectively stable and unstable 
manifolds, that intersect on this critical point. Bra's and 
ket's are forms in the transverse directions and have no 
coordinate dependence in the tangential directions. This 
picture must hold for more intricate flows of hyperbolic 
critical points. The coordinate dependence of the ket of 
the ground state must be related to Lyapunov function 
on the stable manifold. Such function must always exist 
on the stable manifold. 

The closest concept in dynamical systems theory we 
managed to find in the Literature are discussed in Sees. 
3.3-4. of Ref . |32| ■ There, the distributions that localize 
to unstable manifolds of critical points were identified as 
Gel'fand-Schwartz distributions. 

The analysis can be generalized to cases when critical 
points are not isolated and form higher-dimensional crit- 
ical manifolds, M c (Morse-Bott-type case). The analysis 
near each point of M c is split into the directions tangent 
and transverse to M c . In tangent directions, the Hamil- 
tonian is the Laplacian and ground states are harmonic 
forms, |wj}), AmJ^i)) = 0. In transverse directions, if 
conditions are right, the ground state is the one discussed 
previously in the context of an isolated critical point, i.e., 
the ket, |/i}}, is a volume form on stable manifold and is 
localized to the unstable manifold. Provided that |/u))'s 
at different points of M c can be glued together over the 
entire M c , the ground states are \i,M c }) = \cji)) A |/z)). 

Similar situation occurs for perturbative ground states 
around invariant manifolds, Mi. Again, the analysis is 




M,«M 



4/l u JM.» i a » 




FIG. 3: (a) In deterministic limit and under some general 
conditions on the flow, bra's, ((M»|, and ket's, \Mi)), of per- 
turbative ground states around an invariant manifold, Mi, 
represent respectively local stable, M", and unstable, M", 
manifolds intersecting on Mi. (b) An example of perturba- 
tive ground states for a potential-like flow. A stable critical 
point, a, a saddle, 6, and an unstable critical point, c, are 
under consideration. |a)) (black circle) is a maximal-degree 
form that is localized to the critical point, while {{a\ (area 
with horizontal filling) is a constant form with no ghosts and 
represents the entire basin of attraction of this critical point. 
\b)) (think curve aba') and {{b\ (think curve c'bc) are forms 
representing respectively unstable and stable manifolds of b. 
They are forms in transverse directions. |c)) (area with verti- 
cal filling) is a constant form with no ghosts representing the 
basin of repulsion of c, while {{c\ (hollow circle) is a maximal 
degree form localized at this critical point. 



split into tangent and transverse directions. In transverse 
directions, the ground state is the same as in Bott-Morsc- 
type case, |/k)). In tangent directions, one can chose such 
metric go from Eq. (T50)) (if it exists) that A is a Killing 
vector field on Mj. Then, the tangential ground states are 
the invariant harmonic forms also satisfying Ca^i)) = 0. 
The ground states again are as in the Bott-Morse case, 
\i,Mi)) = \uJi)) A |/i)). A toy model of the simplest in- 
variant manifold, a limit cycle, could be a topological 
quantum mechanics on a circle with constant A consid- 
ered in Sec I VI B T1 below. In this case, the ground states 
are from the cohomology of the circle, \uji)) = l,x, i.e., 

|*))o 0) and \^))o } in notations of Scc lVIBll 

For the perturbative analysis, only the combina- 
tions/intersections of bra's and ket's are important. 
These combination are localized to invariant manifolds. 
However, bra's and ket's separately are not localized. In 
fact, perturbative analysis does not say how far bra and 
ket stretch away from the invariant manifold. This pro- 
vides some sort of a freedom in the definition of per- 
turbative ground states. One can use this freedom to 
incorporate the effect of instantons - classical solutions 
that connect different invariant manifolds. The point is 
that instantons provide natural boundaries to stable and 
unstable manifolds of a given invariant manifold. An ex- 
ample of such a definition of perturbative ground states 
is given in FigJSjD, where a potential-like flow in a two di- 
mensional phase space is presented. There, the perturba- 
tive ground states are given for three critical points with 
indices zero (stable critical point, a), one (saddle, b), and 
two (unstable critical point, c). \a)) is a maximal-degree 
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form that is localized to the critical point, while ({a\ is a 
zero-ghost constant form that represents the entire basin 
of attraction of a. \b)) and ((b\ are forms localized to 
respectively the unstable and stable manifolds of b and 
they both have ghosts in the transverse directions. |c)) 
is a constant form representing the "repulsion" basin of 
c, while ((c\ is a maximal degree form localized to this 
critical point. 

As is seen, ket's of perturbative ground states represent 
the instantonic CW-complex of the phase space, while 
bra's represent the dual anti-instantonic CW-complex. It 
is also seen that Q operator acts on perturbative ground 
states just as a boundary operator would have acted on 
the corresponding CW-complex. For example: 



or, similarly, 



Q\b)) = |o» - |a'», 



((&|Q = ((c|-«c'|. 



(60) 



(61) 



In case of a potential flow with isolated critical points, 
perturbative ground states and equalities like the one 
in Eq. (|6"0"|) form Morse- Witten complex, cohomology of 
which must be equivalent to the homology of the instan- 
tonic CW-complex as follows from above discussion. Un- 
der certain conditions, methodology of the Morse-Bott 
complex can be generalized to non-potential flows. The 
proposed view on the perturbative ground states can be 
useful for such generalization. 

From the global point of view, the perturbative ground 
states are not Q-closed as is seen, e.g., from Eq. (|60| . 
This, however, plays no role on the perturbative level 
since, e.g., ((b\Q\b)) — 0. In order to get the global 
ground states beyond the perturbative analysis, one must 
glue perturbative ground states to form globally defined 
stable and unstable manifolds, i.e., consider superposi- 
tion of perturbative ground states that are globally Q- 
closed. 

In the dynamical systems theory, there are theorems 
that establish the conditions, under which the global 
topologically well-defined stable and unstable manifolds 
criss-crossing on invariant manifolds exist. It is natural 
to expect that these conditions are also necessary con- 
ditions for unbroken Q-symmetry in the deterministic 
limit. Those may not be sufficient conditions, however. 
Establishing the conditions of unbroken Q-symmetry in 
the deterministic limit in terms of the properties of the 
flow may turn out to be a non-trivial problem that we 
leave open. We would only like to mention here, that to 
our temporary and mostly intuitive understanding these 
conditions may have a lot to do with the integrability 
conditions. In terms of the spectrum of Fokker-Plank 
Hamiltonian, however, the necessary and sufficient con- 
ditions for unbroken Q-symmetry look fairly clear (see 
FigEb). 

As a closing part of this section, we would like to point 
out the equivalence of Eq.(60) and the supersymmctry 
operator in stochastic networks discussed, for example, 



in Ref. 28]. There, the supersymmetry relates probabil- 
ity currents and probability densities defined on a net- 
work. The direct analogy with our case is as follows: the 
supersymmetry is the exterior derivative, the network is 
the instantonic CW-complex, and, say, \b)) is the current 
between \a)) and \a')). 



V. UNBROKEN Q-SYMMETRY 
A. Langevin SDEs 

The important class of models with unbroken Q- 
symmetry is the Langevin SDEs. They represent purely 
dissipative dynamics. This is the most studied class of 
dynamical systems in the context of W-TFTs. The mod- 
els are quasi-Hcrmitian and have_a Hermitian represen- 
tation known as Witten models. [9| As we discuss in this 
subsection, Q-symmetry is never broken for Langevin 
SDE's (for equilibrium dynamics) even in the presence 
of noise. 

Langevin SDEs are those with a potential flow: 



Ai(x) = g ij 6 vi(x) W, 



(62) 



where W is some functional called superpotential. All 
the eigevalues of the corresponding Fokker-Plank Hamil- 
tonian are real. The Hamiltonian can be brought to a 
Hermitian form by the similarity transformation 



H^H L = Hl = fil /2 Hfi L 1/2 , 



(63a) 



where t)l = s 2W is the Hilbert space metric discussed 
in Sec lIII PI The corresponding transformation of the 
Hilbert space and the supercharges is 



«*| -► <*| L =«tt|C 



-1/2 



l*» -► I*)l = ^ /2 |*» 



(63b) 
(63c) 

Q -► Ql = Q+ I X l (x){5 v , (x) W), (63d) 

J X 

Q -► 4 = O f - fs xHx) g ij (S vHx) W). (63e) 

J X 

It is easy to see in this hermitian form that the model has 

yet another nilpotent supersymmetry charge, Ql,Ql = 

0, [H, Qi]- — 0, which is the Hermitian conjugate and 

the time-reversal companion of Ql = 0L- The model 
became N=2 supersymmetric with the Hamiltonian 

2H L = [Q L ,Ql} + = Ql = Ql (64) 

where the two Hermitian supercharges are Q\ = Ql + Ql 

and Q\ = z(Ql - Ql)- 

In this representation, all the eigenvalues are real and 
non-negative, £ n — T n : 



2r n = |Q L |n) L | 2 + |Q L |n) L | 2 >0. 



(65) 



13 



For ground states we have r„ = so that: 

Qi»l = 0,4» l = 0. 



(66) 



Now it follows that Q-symmetry, together with the N=2 
supersymmetry, can not be broken if at least one physical 
ground state exist. Such ground state does always exist. 
It has the from of Eq. (|27|) : 



(67) 



|0)) = (*Po),«OH(*|0)})*77l = 1. 



Here we turned back to the original bi-orthogonal basis 
in the Hilbert space in order to emphasize the physi- 
cal meaning of this ground state - the stationary prob- 



ability distribution given by Pq 



-2W 



We only con- 



sider physically meaningful superpotentials, W. In par- 
ticular, W must be bounded from below. As is seen, 
H cnv Po = 0, with the conventional Fokker-Plank op- 
erator from Eq. ([29|) . or, equivalently, H\0)) = 0. Un- 
conditional existence of this ground state suggests that 
Q-symmetry is never broken for Langevin SDE's with 
physically meaningful W's. 

The above discussion suggests that a model must have 
a non-potential part in its flow for the Q-symmetry 
to be spontaneously broken. In physical terms, to 
break Q-symmetry of a Langevin SDE one has to ex- 
ternally "drive" the system through its phase space by 
a non-potential flow. This driving is actually a well- 
documented condition for SOC dynamics. [11| 

It is possible that there exist other classes of models 
with unbroken Q-symmetry even in the presence of noise. 
One of the possible candidates is the conservative models. 
If it turns out to be true, one would say that a nearly con- 
servative model must possess dissipative part in its drift 
term for Q-symmetry to be broken. This is actually one 
of the conditions for the observation of chaotic behavior 
in nearly conservative models. Therefore, it is reasonable 
to expect that in some subclass of stochastic conservative 
models Q-symmetry is also unbroken. 



B. Quantum excitations vs. stochastic fluctuations 

Here we would like to address the following subtle 
point. It is known that W-TFT's with unbroken Q- 
symmetry do not possess physical excitations (see, e.g., 
Refs. 0, Q as well as Sec lIIIF[) . This may sound sus- 
picious: how can a theory with no quantum excita- 
tions represent a stochastically fluctuating system such 
as a Langevin SDE? To get around this controversy, in 
Ref.[29( it was wrongfully conjectured that the PSW 
stochastic quantization is applicable only to stochastic 
systems with slow colored noises that do not provide high 
enough frequencies for the system to fluctuate. This is 
not true - PSW method is applicable to all dynamical 
models. 

The resolution of the seeming controversy is this. Con- 
sider again Langevin SDEs as an example. They have 
only one physical ground state of trivial ghost content 



- the stationary probability distribution from Eg. (1671) . 
This, however, does not mean that there are no stochas- 
tic fluctuations. The very necessity to use the probability 
distribution for the description of the model suggests that 
the bosonic variables fluctuate. These fluctuations, how- 
ever, are in the sense of the zero quantum fluctuations 
within a ground state. 

Furthermore, W-TFTs with unbroken Q-symmetry 
provide perfectly defined correlators that represent the 
stochastic fluctuations in the bosonic variables. For 
example, the one-loop propagator for the fluctuating 
bosonic fields near an isolated stable critical point, 
ifo, of a (0+1) theory (not necessarily Langevin) is 

«^ |V(*)V(0)ko)) = {&D)~\ where D = 5)d t + 
d ip jA l (ipo) and \<po)) 1S the one-loop ground state at 
ipo. These correlators are not BPS observables and thus 
bclongto the general class of observables studied in 
Refs. flO| . The above propagator is exactly the one- loop 
propagator that would follow from the more wide-spread 
MSR approximation to stochastic quantization [higher- 
order perturbative corrections will be different, however, 
due to the neglect of the contributions from the virtual 
ghosts in MSR]. 



VI. SPONTANEOUSLY BROKEN 
Q-SYMMETRY 

In thinking about models with spontaneous Q- 
symmetry breaking, it is an important question on which 
level the symmetry is broken. For any global continu- 
ous symmetry, there are three such levels - the mean- 
filed level, the perturbative level also known as quan- 
tum anomaly, and by the dynamical condensation of 
instantons and anti-instantons. The last possibility is 
particularly important for supersymmetries such as Q- 
symmetry. [6|] The reason for this is that certain class 
of supersymmetries, including Q-symmetry, can not be 
broken on the perturbative level. This is related to the 
concepts of the cancelation of bosonic and fermionic de- 
terminants, of the localization of path-intcrgals to clas- 
sical solutions in W-TFT's (see, e.g., Ref. [J]), and to 
non-renormalization theorems. 

This leaves out only two possibilities: Q-symmetry can 
be broken either on the mean-field level (which in our 
case corresponds to the deterministic limit) or by the dy- 
namical condensation of instantons and anti-instantons. 
Below we argue that these two situations represent re- 
spectively chaotic behavior (Sec. IVIBp and Intermit- 
tent/SOC dynamics (Sec. IVTCj) . 

Dynamical (anti-)instanton-mediated breaking of Q- 
symmetry must not exist in the deterministic limit be- 
cause of the disappearance of anti-instantons. Never- 
theless, in non-equilibrium situations such as quenches 
instantons are not required to be matched by anti- 
instantons and consequently instanton(s) alone can break 
Q-symmetry even in the deterministic limit. [lOj This pos- 
sibility is briefly discussed in Sec I VI Dl Let us first, how- 
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ever, address some general aspects of Q-symmetry broken 
dynamics. 



A. Qualitative difference of Q-broken dynamics 

Models with unbroken and spontaneously broken Q- 
symmetry must have qualitatively very different dynam- 
ics. This difference has a few separate points. Allow us 
address some of them here. 

First, spontaneously broken Q-symmetry results in the 
liberation of quantum dynamics in the following sense. If 
one randomly chooses a wavefunction, |\&)), and propa- 
gates it in time long enough, he will have: 



!*»(*) 



-iE„t 



|n))6W 



n»«n|tf)). 



(68) 



Now it is seen that in case of unbroken Q-symmetry, 
when the only physical states are the ground states with 
E n = (see FigO)), the quantum evolution in Eq. (j6"8")) 
stops. In contrary, if Q-symmetry is spontaneously bro- 
ken, the time evolution will never stop. 

If one would want to construct an effective held theory 
out of a Q-broken model, he could as well try to write 
down a theory that covers the physical states only. Such 
theory would necessary be Hcrmitian. In other words, 
on physical states the Fokker-Plank equation becomes a 
Schodinger equation without Wick rotation, i.e., preserv- 
ing the original meaning of time. This line of thinking 
can be useful for establishing a long-suspected connec- 
tion between chaotic and quantum dynamics, (see, e.g., 
Rcf. 30]) W-TFT's may provide a firm basis for such 
connection. 

Second, Goldstone theorem ensures that if a global 
continuous symmetry is spontaneously broken then there 
must exist a gapless particle. In case of Q-symmetry, the 
gapless particles are ghosts called Goldstinos. Goldstinos 
provide the system with long-range spatiotemporal corre- 
lations in some observables related to them. These long- 
range correlations are known under various names such as 
sensitivity to initial conditions, self-similarity, algebraic 
statistics of avalanches in SOC dynamics, 1/f noise, non- 
Markovian memory etc. Such long-range correlations 
must not exist in models with unbroken Q-symmetry that 
are thus could be called Markovian dynamics. 

Third, the time-reversal symmetry must also be bro- 
ken in Q-broken phases as was discussed in the end of 
Scc lIIIFl This must be related to the concept of irre- 
versibility of, e.g., chaotic dynamics. 



B. Chaotic dynamics 

As it follows from the previous subsection, models with 
spontaneously broken Q-symmetry must have qualita- 
tively very different behavior. In deterministic limit, 
there is only one class of dynamics general enough to 




FIG. 4: (a) The spectrum of a toy model considered in 
Sec lVIBTl In the deterministic limit, e — > 0, the Fokker- 
Plank eigenvalues tend to lay on the imaginary axis, as though 
Q-symmetry is broken. This situation must also occur in 
other conservative models. As is explained in the text, this sit- 
uation must not be viewed as the mean-field level Q-symmetry 
breaking, (b) The generic phase diagram of a stochastic dy- 
namical system. It consists of three phases. The first is the 
chaotic phase, where the topological symmetry is broken on 
the mean-field level. The second phase is intermittent chaos 
or self-organized criticality, where the topological symmetry 
is dynamically broken by the condensation of instantons and 
antiinstantons. SOC is not a "critical state" but rather is a 
full-dimensional phase. However, in the deterministic limit 
when antiinstantons disappear, the SOC phase collapses into 
the critical "edge of chaos". Chaotic and SOC phases must 
possess spatiotemporal self-similarity due to the existence of 
gapless goldstinos. The temporal aspect of the self-similarity 
represents dynamical non-Markovian long-term memory. By 
the same token, the phase with unbroken topological symme- 
try can be called Markovian dynamics (MD). 



be associated with the mean-field level spontaneous Q- 
symmetry breaking on the W-TFT side. This class is the 
deterministic chaos. 

In the deterministic limit, the effects of quantum- 
mechanical tunneling (anti-instantons) between pertur- 
bative ground states can be neglected. Perturbative 
ground states discussed in Sec HVI become a good approx- 
imation for unbroken Q-symmetry cases. Then, the rhs 
of Eq. ([rj5|) for a randomly chosen initial probability dis- 
tribution (a wave- function of the maximal degree, \^))) 
is the sum over all the stable attracting invariant mani- 
folds, while the overlapping coefficients, {{n\^)), are the 
integrals of |\1/)) over the corresponding attracting basins. 
The absence of quantum evolution means that the deter- 
ministic model with unbroken Q-symmetry will always 
find its attractor and stay there forever. In contrary, 
when Q-symmetry is broken and the evolution nether 
stops, it must mean that a model is not capable of finding 
its topologically well-defined attractor - a feature perti- 
nent to chaotic behavior. 

Yet another way to see that mean-field Q-symmetry 
breaking must be associated with chaos it through con- 
sidering models with fractal strange attractors. In such 
cases, one will certainly find it difficult to represent such 
an invariant "manifold" (formally a fractal is not a topo- 
logical manifold) as an intersection of forms, that the 
ground states are supposed to be in the unbroken Q- 
symmetry models. The first problem here is the differen- 
tiability issue. The second problem is more fundamental. 
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The point is that states can only be forms of integer de- 
gree. Such forms can only represent integer dimensional 
manifolds but not fractals. 

Looking at deterministic chaotic models as at the 
mean-field level Q-symmetry breaking suggests why they 
must exhibit long-range correlations in observables re- 
lated to ghosts. These correlations is the W-TFT way to 
explain the high sensitivity to initial conditions in chaotic 
models. In dynamical systems theory, this sensitivity is 
often considered as one of the necessary conditions for 
the system to be chaotic. On the W-TFT side, the sensi- 
tivity is the consequence of the Goldstone theorem. Fur- 
thermore, time-reversal symmetry of deterministic chaos 
must also be broken. This is most probably the W-TFT 
way of encoding the temporal irreversibility of chaos. [331 ] 

The above reasons are only indications on the equiv- 
alence between deterministic chaos and mean-field level 
spontaneous Q-symmetry breaking. A rigorous prove of 
this equivalence has obstacles. One of them is the absence 
of the universal definition of deterministic chaos. From 
this perspective, one can as well think that W-TFT's ap- 
proach provides its own version of this definition. 

If Q-symmetry is broken in the deterministic limit, i.e., 
on the mean-field level, the addition of noise will not re- 
store it. In fact, as is discussed in the next subsection, 
the noise can only break Q-symmetry if it was not bro- 
ken on the mean-field level. Stochastic models with Q- 
symmctry broken on the mean-field level can be called 
stochastic chaos. 



1. Refined definition of deterministic chaos 

So far, the W-TFT definition of chaos is a model 
with Q-symmetry spontaneously broken on the mean- 
field level. This definition needs further clarification. The 
reason for this is seen from the following example. 

Consider a dynamical model on a circle, S 1 , of circum- 
ference, L, with constant flow, A. The operators defining 
the Fokker-Plank Hamiltonian through Eq. (l26aj) are: 



Q = xd< P ,Q = d x (-ed (p -2A). 



(69) 



The eigenstates and the spectrum are known exactly: 
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££' = ££> = ek z /2 - 2ikA, 



(70) 
(71) 



where ip £ [0, L] and the momentum is kL/(2ir) 6 Z. As 
is seen in FigfJJ., the spectrum looks like that of broken 
Q-symmetry in the strict deterministic limit, e = 0. On 
the other hand, for any small but finite noise, < e <C 
1, the spectrum is that of unbroken Q-symmetry. Such 
situation must appear in all conservative models. 

To get around this seeming controversy, we have to 
refine the definition of the mean-field level Q-symmetry 
breaking and thus that of chaos. Chaotic system is such 
that its Q-symmetry is spontaneously broken in the limit 



of infinitely weak but finite noises. This weak but finite 
noise is related to the concept of sensitivity to initial con- 
dition on the dynamical systems side. To see if a deter- 
ministic model is sensitive to initial conditions one must 
consider infinitely close but still different initial points in 
the phase space. This can be interpreted as the addition 
of infinitely weak but still finite uncertainty (noise). 



C. Intermittency/SOC 

In stochastic models, there is yet another mechanism 
for spontaneous Q-symmetry breaking. This mechanism 
is the dynamical condensation of instantons and anti- 
instantons. 0, @, 0] We have already encountered in- 
stantons in SeclIVI in the discussion of the ground states 
in the deterministic limit. Let us, however, briefly rec- 
ollect again on what these tunneling processes are. In- 
statnons correspond to classical solutions of SDE (solu- 
tion of a corresponding DDE) that start on one invariant 
manifold and end at another. Thus, instantons always 
lead from " less stable" invariant manifolds to more stable 
ones. Anti-instantons, in turn, are time-reversed instan- 
tons. They lead from more stable invariant manifolds to 
less stable ones. Anti-instantons are essentially the mo- 
tion against the flow, which is only possible in the pres- 
ence of noise. Matrix elements of anti-instantons always 
come with exponentially small factors that disappear in 
the deterministic limit. A hydrodynamical example of 
an instanton/antiinstanton is the processes of annihila- 
tion/creation of (pairs of) vortices or vortex lines. 

(Anti-)instantons intrinsically break Q-symmetry as, 
in particular, is seen in Eg. (1601) . from which it follows 
that the matrix element of the expectation value of a Q- 
symmetry operator on the (6a)-instanton (see FigjS})) is 
non-zero: ((a\Q\b)) = 1. 

In equilibrium situations, condensation of instantons 
must be accompanied by the condensation of antiinstan- 
tons. Indeed, for the instantonic processes of annihilation 
of vortices to be happening forever they must be com- 
plemented by the antiinstantonic processes of creation of 
vortices. In other words, separately each (anti-)instanton 
is ultimately non-equilibrium processes leading from one 
state to another. Therefore, in equilibrium situations, 
only configurations of "equal number" of instantons and 
anti-instantons that lead from a state to itself exist. 
Thus, in equilibrium, (anti-)instanton-mediated dynam- 
ical Q-symmetry breaking can only happen in stochastic 
but not deterministic models. 

The dynamics in a phase with the (anti-)instanton- 
mediated Q-symmetry breaking must look as an infinite 
series of jumps (e.g., avalanches of sandpile models) be- 
tween different invariant manifolds. If the invariant man- 
ifolds are points in the phase space (solitonic configura- 
tions, patterns etc.), we can talk about the dynamics as of 
a sequence of sudden changes in solitonic configuration. 
This is nothing else but the physical picture of SOC. [ll| 

In a more general case, these jumps happen not be- 
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tween points but between more intricate invariant mani- 
folds (e.g., limit cycles). In this case, the term SOC must 
be substituted by a more general concept of Intermit- 
tency (see. e.g., Ref . |34| ) - in dynamical systems theory a 
system is called intermittent if it is subject to infrequent 
variations of large amplitude. These variations ((anti- 
)instantons) separate periods of different (quasi-)periodic 
behaviors. Roughly speaking, the difference between In- 
termittency and SOC is that the jumps happen between 
limits cycles or more complicated higher-dimensional sta- 
ble sets for the former and between stable points for the 
later. Thus, SOC is in a sense a subclass of intermit- 
tent behavior. We do not see, however, more fundamen- 
tal difference between SOC and Intermittency and thus 
sometimes use these terms interchangeably. 

The conclusion we just arrived is that (anti-)instanton- 
mediated Q-symmetry broken phase must be associated 
with Intermitency/SOC. [29| On the phase diagram, in- 
tcrmittcnt/SOC phase must occupy a region between 
chaotic and Markovian phases. Importantly, Intermit- 
tency/SOC is not a critical state. It is a rightful full- 
dimensional phase, just as chaos and Markovian systems. 
This explain why by moderate variation of parameters 
of the model one can not destroy intermittent/SOC be- 
havior. In the literature of SOC, this fact is sometimes 
attributed to the mysterious tendency of SOC systems to 
" evolve" to a critical " edge of chaos" . W-TFT picture of 
Intermittency/SOC resolves this issue. 

Only in the deterministic limit, anti-instantons vanish 
and intermittent/SOC phase collapses into the "edge of 
chaos" (see, e.g., Ref. |35j ) . The above discussion is sum- 
marized in Fig.((3)D) where chaotic, intermittent/SOC, 
and Markovian phases are schematically presented. 

The proposed phase diagram complies with previous 
studies. One of the examples is the phase diagram of wa- 
ter moving between two rotating cylinders. [36j The Cou- 
ette (laminar) flow phase can be identified with Marko- 
vian phase, the featureless turbulence with chaos, while 
the intermediate region, where the dynamics is domi- 
nated by creation/annihilation of solitons in the form 
of vortices and vortex lines, must be identified with in- 
termittency/SOC. Yet another example is neuroscience, 
where the three phases in Fig.((3t>) are often called re- 
spectively supercritical, critical, and subcritical phases. 



D. Non-equilibrium dynamics 

Let us also briefly address non-equilibrium situations. 
In these situations, the temporal boundary conditions 
are not periodic and the system is allowed to end up in 
a state different from the initial state. Physically, this 
corresponds to quenches for instance. Quench dynam- 
ics can be assumed deterministic and it can be looked 
upon as a composite instanton leading from one of the 
metastable configurations to one of the stable configu- 
rations. Barkhausen-like effects (see, e.g., Ref.|38| and 



Refs. therein) such as crumpling paper (see, e.g., Ref.|39j 
and Refs. therein) can effectively be viewed as "slow" 
quenches where the time evolution is due to varying ex- 
ternal parameter. Gradual variation of the external pa- 
rameter (e.g., magnetic field in Barkhausen effect) makes 
previously stable configurations unstable ones thus ini- 
tiating instantons, that look like a sequence of sudden 
jumps. As follows from the next paragraph, these mod- 
els must also exhibit "power-laws" which is the reason 
why they also sometimes attributed to the SOC family. 
Mathematical aspects of non-equilibrium deterministic 
dynamics in potential flows (Langevin DDE's) were stud- 
ied, e.g., in Ref. [lfj| . There, it was found that the cor- 
responding theory is a log-conformal theory. Roughly, 
one can attribute the long-range log-conformal correla- 
tions to the Q-symmetry spontaneously broken by a given 
composite instanton. Interestingly, a quench must not 
necessary be across some-other-symmetry breaking phase 
transition in order to exhibit long-range correlations, e.g., 
Barkhauscn-likc dynamics. It is also interesting that for 
quenches one can straightforwardly use BPS observables 
that may provide certain topological invariants in topo- 
logically non-trivial theories. It is tempting to believe 
that so-calculated topological invariants may have some- 
thing to do with the probabilities of the system to end 
up in one or the other competing stable configurations. 



VII. CONCLUSION 

In this paper, it is shown that the most general 
stochastic quantization procedure applied to any stochas- 
tic or deterministic continuous-time (partial) differential 
equations leads to a Witten-type topological field the- 
ory - a theory with global topological supersymmetry. 
This topological symmetry must be perturbatively sta- 
ble and consequently can be spontaneously broken only 
either on the mean-field level by, say, fractal invariant 
sets, or dynamically by the condensation of instantons 
and anti-instantons. According to this, we propose a 
generic phase diagram given in Fig. ([!]}). The mean- 
field level spontaneous Q-breaking must be associated 
with chaotic behavior. The instanton-anti-instanton me- 
diated breaking of topological symmetry corresponds to 
intcrmittcncy/SOC - a full dimensional phase surround- 
ing chaos. Its full-dimensionality explains why algebraic 
correlations of SOC dynamics are robust against mod- 
erate variation of the parameters of a system. Full- 
dimensional intermittency/SOC phase exists, however, 
only in stochastic models. In deterministic models, when 
anti-instantons disappear, intermittency/SOC phase de- 
generates into the critical " edge of chaos" between deter- 
ministic chaos and deterministic Markovian dynamics. 

According to Goldstone theorem, both phases with 
spontaneously broken topological symmetry must exhibit 
long-range correlations. These correlations are related 
to such well-established concepts as self-similarity, sensi- 
tivity to initial conditions and/or non-Markovian scale- 
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free memory, power-law statistics of avalanches, algebraic 
power-spectrum in turbulence etc. Accordingly, dynam- 
ical systems with unbroken topological symmetry can be 
called Markovian in that sense that they do not exhibit 
long-range dynamical memory so one can always look at 
the system at such temporal scale that the dynamics will 
look effectively Markovian. [40] 

The phase diagram and the interpretation of the phys- 
ical essence of the three phases is the only result of this 
paper. There is certainly much more in the W-TFT of 
dynamical systems to be understood in the future. In 
particular, in this paper we never used one of the most 
fascinating possibilities that W-TFT's offer - the possibil- 
ity to calculate certain topological invariants as expecta- 
tion values of BPS observables on instantons. One of the 
candidate ways to exploit these observables (probably in 
combination with their anti-instantonic counterparts) is 
in the form of order parameters for Q-brokcn phases. In 



fact, a rigorous way for construction of effective field the- 
ories is one of the important advancements that W-TFT 
may lead to. W-TFT of dynamical systems may also hold 
the key to some novel topological aspects of probability 
theory. 
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